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d = -If t 1 " cos ( 2lt ^ T )3 (H) 

where the constant of integration has been evaluated by setting d = 0 
at T =: 0. 


The position of the wake center at x for successive instants of 
the time t may he derived from consideration of the following displace- 
ment sketch: 



A line element released from the line source at time ^ q reaches the 
axial station x after a residence time x/u* The lateral displacement 
y from the flow-axis plane is given by 

y = d(x) - d(x 0 ) = d(x) - d(x - (12) 


where d(x) is d evaluated at time x and d(x - x/u) is d at time 
x - (x/u). By combining equation (ll) with equation (12), the time 
history of the position of the center of the molecular-diffusion wake at 
any axial station x is 


y = 



(2it/x)J - jl 


cos 


^2jt/x - 2 ji/ 



(13) 
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5 By William R. Mickelsen and Lionel V. Baldwin 


SUMMARY 

The processes of fuel-air preparation, combustion, and exhaust-gas 
mixing in jet -engine combustors depend to a large degree on aerodynamic 
mixing. Since combustors commonly have intense sound fields, it is of 
interest to investigate the effect of sound on the fundamental mixing 
process. This report describes a theoretical and experimental investi- 
gation of the aerodynamic mixing by a standing sound wave downstream 
from a continuous line source of heat. 

By a kinematic analysis of the motion of the molecular -diffusion 
wake, equations are derived for the time variation of temperature and the 
time-mean temperature at points throughout the mixing region. The 
analysis shows that standing sound waves displace the diffusion wake in 
a manner similar to the displacements of a flag waving in a harmonic 
mode. The diffusion-wake displacement has nodal points downstream from 
the line source at distances which are integer multiples of the ratio 
(stream velocity)/ (sound-wave frequency). The maximum wake displacement 
occurs at the antinodal distances and is equal to x rms transverse 

velocity fluctuation)/(rt x sound-wave frequency). The theoretical 
analysis considers two general cases; (a) small and (b) large geometrical 
deformations of the diffusion wake by the transverse sound waves. The 
analysis shows -that the deformations are characterized by stretching of 
the wake and by the relative magnitudes of the wake radius of curvature 
and thickness. The stretching is greater for high sound intensities. 

The curvature is greater for high sound intensities and high sound 
frequencies. When the wake suffers negligible deformation, the sound 
field contributes to the mixing in the time-mean sense; but the instan- 
taneous, or local spatial, structure of the temperature field is unaf- 
fected. For the case of appreciable wake deformation, the temperature 
field is substantially changed in both the time-mean and instantaneous 
senses. 
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Experimental measurements confirmed the theoretical analysis for the 
case of negligible deformation of the diffusion wake, but no measurements 
were made for the case of appreciable deformation. The experiment was 
carried out in a low-turbulence airstream with superimposed transverse 
velocity fluctuations associated with a standing sound wave. The mixing 
region downstream from the line source of heat was at a pressure node of 
the standing sound wave so that pressure fluctuations were negligible. 

The amplitude of the transverse velocity fluctuations corresponded to a 
sound pressure level of 147 decibels, and the wave frequency was 104 
cycles per second. Time -mean temperatures were measured with a thermopile 
probe. Instantaneous temperatures were measured by two different methods: 
resistance thermometry and hot-wire anemometry. The velocity fluctuations 
in the sound field were measured by hot-wire anemometry. 

In order to facilitate application of the results to practical 
problems, a special section is included in which equations and methods 
are summarized. 


INTRODUCTION 

Aerodynamic mixing is an important factor in fuel-air preparation, 
combustion, and exhaust-gas mixing in jet-engine combustors. Because of 
its importance, it is of interest to investigate aerodynamic mixing under 
flow conditions commonly found in combustors. Although jet -engine com- 
bustors often have extremely intense sound fields, little attention has 
been given to the effect of sound-wave disturbances on aerodynamic mixing. 
The research described in this report covers a fundamental investigation 
of aerodynamic mixing in a low-turbulence airstream having a s ingle - 
frequency, plane, transverse, standing sound wave. 

Mixing in fluid streams proceeds locally by molecular motion 
(molecular diffusion) and on a larger scale by the turbulent velocity 
fluctuations, which transport fluid volumes composed of many molecules. 

It is reasonable to expect that the velocity fluctuations arising from 
the passage of sound pressure waves contribute to the mixing process in 
a manner somewhat similar to that of turbulence. The data in reference 1 
show that the approach stream in a smooth-burning ram-jet combustor con- 
tained sound and turbulence velocity fluctuations of the same order of 
magnitude; the data of reference 2 show that screeching combustors have 
velocity fluctuations due to sound that may be several orders of magnitude 
larger than the turbulent velocity fluctuations. 

Since the velocity fluctuations responsible for aerodynamic mixing 
are usually slow compared with the flow velocity through the combustor, 
the concentration of heat or mass in mixing regions fluctuates with time 
and varies irregularly from point to point in space. In turbulent mixing, 
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concentration fluctuations in time and space may also be due to the finite 
range of turbulent eddy sizes, since the smallest turbulent eddies are of 
a much larger scale than molecular mean free path. Experiments in turbu- 
lent mixing (refs. 3 and 4) have shown that concentration fluctuations 
can be of nearly as great a magnitude as the mean concentration at points 
throughout the mixing region. For these reasons, it is evident that 
mixing should be considered in the instantaneous sense as well as in the 
time-mean sense. Mixing in the instantaneous sense can be described by 
the time history of concentration fluctuations at a point in space, and 
mixing in the time-mean sense can be described by the long-time average 
of concentration at a point in space. 

In order to investigate the fundamentals of the aerodynamic mixing 
due to sound, an experimental and theoretical study was made of the 
mixing downstream from a line source of heat in an intense sound field. 

The turbulence in the airstream was kept at a negligible level so that 
the mixing process was due only to the sound field and molecular diffu- 
sion. To prevent vortex shedding, the electrically heated wire that 
formed the line source had a small diameter. The center plane of the 
two-dimensional mixing region was placed at the velocity antinode and 
pressure node of a plane standing sound wave so that the sound pressure 
fluctuations were negligible and the sound disturbance velocities were 
transverse to the stream direction as shown in the following sketch: 


Transverse velocity fluctuations imposed by standing sound wave 

The magnitude of the transverse velocity fluctuations corresponded to a 
sound field having a sound pressure level of 147 decibels. The trans- 
verse velocity fluctuation was nearly sinusoidal and had a frequency of 
104 cycles per second. The instantaneous temperature was measured by a 
resistance thermometer and a hot-wire anemometer, and the transverse 
velocity fluctuations were measured by hot-wire anemometry. The time- 
mean temperatures were measured with a thermopile probe. 


velocity 


Stream 


Lin 
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By a kinematic description of the mixing region, equations are de- 
rived for the time variation of temperature and for the time-mean tem- 
perature at points throughout the mixing region. This theoretical analysis 
is made for two cases of interest: (a) negligible geometrical deforma- 

tions of the diffusion wake by the sound field and (b) appreciable 
deformations of the diffusion wake. For the case of negligible deforma- 
tion of the diffusion wake, the theoretical analysis is compared with 
measured values of the instantaneous and mean temperatures at various 
points downstream from the line source of heat. By means of the theo- 
retical analysis, a comparison is made between the cases of negligible 
and appreciable wake deformation. 


THEORETICAL ANALYSIS 

The analysis of the temperature field in the mixing region is based 
on a kinematic description of the displacement of a molecular -diffusion 
wake by a transverse, plane, standing sound wave. A pictorial represen- 
tation of the wake displacement is shown in figure l(a). As shown in 
the following analysis, the displacement of the wake is similar to the 
displacement of a flag waving in a harmonic mode. The center of the 
wake has a sinusoidal shape; and the molecular diffusion proceeds as 
the wake elements pass downstream, as shown by the temperature profiles 
in figure l(a). From the kinematic analysis, expressions are derived 
for the time variation of temperature and for the time-mean temperature 
at fixed points throughout the mixing region. 


Molecular -Diffusion Wake 

The molecular diffusion of heat in a fluid stream at constant pres- 
sure is described by the Fourier -Poisson equation (ref. 5) derived from 
the Fourier -Biot heat-flux law. The Fourier -Poisson equation is 

§ - sr f ir) t 1 - 1 < 2 - 3 > ^ 

where Dh/Dt is the total derivative of the volume specific enthalpy 
h, which is defined by 


T 


L am 


h = 


pep dT 


( 2 ) 


4022 


ZZOf 


NACA TN 3760 


5 


and the coefficient of thermal diffusivity a is 

k 



( 3 ) 


(Symbols are defined in appendix A. ) 


For the case of molecular diffusion of mass, an equation similar to 
equation (l) can be derived from Fick's mass-flux law: 



(i = 1,2,3) 


( 4 ) 


where is the molecular-diffusion coefficient and ^ is the concen- 

tration of the diffusing gas. Since equations (l) and (4) are identical 
except for the coefficients a and D^, their solutions for specific 

problems are also identical except for the coefficients. Because of the 
close similarity of equations (l) and (4), the following theoretical 
analysis applies in essence to the mixing of fluids with different heat 
contents or of different physical properties. 

If the coefficient of thermal diffusivity is essentially constant, 
and if a molecular spreading coefficient to is defined by 

'H 

cu = a dt (5) 

Jo 


then equation (l) may be written 


Dh _ d 2 h 
Da> Sx^ 2 


(i = 1,2,3) (6) 


If the fluid field is at rest, or if a coordinate system is used that 
moves along with the stream, then equation (6) can be written in terms 
of a Cartesian coordinate system (?, rj, £) : 


dh _ ^h + c^h + c^h 
^ c £ 2 dr ] 2 S ? 2 


As shown in reference 6, the solution of equation (7) for an instantan- 
eous line source lying along the £-axis and emitted at t = 0 is given 
by 


h 


= iZi 

4rtcu 


e 


g 2 + 11 2 
4u> 


( 8 ) 
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where H/l is the heat per unit length emitted at t = 0. Equation (8) 
is used in the subsequent derivation of an exact model for instantaneous 
temperature in the mixing region. 

By assuming that diffusion in the stream direction is negligible, 
the following approximate equation for diffusion from a continuous line 
source in a stream flowing with velocity U in the x-direction is de- 
rived in reference 7 : 


h = - 9 l> - . e" (9) 

U 

where T] is the lateral distance from the molecular -wake center to a 
fixed point in space, q/Z is the rate of heat release per unit length 
at the line source, and oi^ = aX/U. Equation (9) is used in the subsequent 

derivation of an approximate model for instantaneous temperature in the 
mixing region. The temperature difference AT = T - Tam can be deter- 
mined from the volume specific enthalpy h by means of tables or graphs 
based on equation (2). 

Kinematics of Molecular -Diffusion Wake in Transverse Sound Field 

A kinematic description of the motion of the molecular -diffusion 
wake is necessary to subsequent derivations of the instantaneous and 
time-mean temperatures in the mixing region. This section describes the 
wake motion in terms of the time T associated with the velocity fluc- 
tuations of the standing sound wave and of the distance x downstream 
from the line source of heat. 

The transverse velocity fluctuation associated with the sound wave 
is assumed to have a sinusoidal form given by 

v = sin (2 ji/t) (l0) 

where *y/v2 is the root-mean-square value of the velocity wave form, / 
is the sound-wave frequency, and T is time. 

A fluid line element in the stream has a transverse motion caused 
by the velocity fluctuation v and a resulting transverse displacement 
d. The time history of the displacement d of a line element in the 
center of the molecular -diffusion wake is given by the integral of 
equation (10) : 
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(H) 


where the constant of integration has been evaluated by setting v = 0 
at x = 0. 

The position of the wake center at x for successive instants of 
the time x may be derived from consideration of the following displace- 
ment sketch : 


A line element released from the line source at time x reaches the 
axial station x after a residence time x/u. The lateral displacement 
y from the flow-axis plane is given by 


where d(x) is d evaluated at time t and d(x + x/u) is d at time 
x + (x/u). By combining equation (ll) with equation (12), the time 
Bi story of the position of the center of the molecular-diffusion wake at 
any axial station x is 



x 



( 12 ) 


V ■ {J 1 - 003 ( 2 *A + 2"/§)l - [1 - (2*A)]} (13) 
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which reduces to 


y , JLp. sin (^) sin (*„ + #) 


( 14 ) 


Figure l(h) shows the position of the molecular-diffusion-wake 
center at various times through one cycle of transverse velocity fluc- 
tuation. The wake positions are shown only through two nodal points, 
hut theoretically they continue on downstream in the same cyclic fashion. 
An interesting feature of figure l(b) is the existence of definite nodal 
points at /x/u =1,2, . . . The reason for these nodal points is shown 
by the sketch in the preceding paragraph. A fluid line element leaving 
the line source at any time t has zero lateral displacement from the 
flow axis if its residence time x/u is any multiple of l//. Also of 
interest is the maximum value of the wake center displacement at the 
antinodal planes /x/u =0.5, 1.5, . . . This maximum value is equal to 


center are dependent only on the intensity and frequency of the trans- 
verse sound wave. 


The description of the instantaneous temperature at points in the 
mixing region is of interest both as an important factor in aerodynamic 
mixing and as a preliminary step in the evaluation of the time-mean 
temperature. A physical picture of the instantaneous temperature may be 
obtained by a study of the molecular -wake -position diagram shown in 
figure l(b ) . As the molecular wake approaches a fixed point, say P]_, 

the temperature at the point rises because of the temperature profile of 
the wake as illustrated in figure l(a). When the wake center reaches 
the point, the temperature is at the peak value of the wake profile. 

Then, as the w^ke passes on away from the point, the temperature decreases 
to a minimum which is reached when the wake is at its farthest position 
from the point. Because the lateral motion of the wake is cyclic, the 
time history of the instantaneous temperature will also be cyclic with a 
period equal to that of the lateral velocity fluctuations associated with 
the standing sound wave. 

Typical time histories of instantaneous temperatures are shown in 
figure 2. The three time histories represent the temperature fluctua- 
tions at points P-j_, P 2 > and P 3 in figure l(b), and were calculated 

from a theoretical equation derived later in this section. Consideration 
of figures l(a) and (b) shows that the temperature at the nodal lines 
will be steady in time. Temperatures at other points in the nodal planes 



that the maximum lateral excursions of the wake 


Instantaneous Temperature in Mixing Region 
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(for instance, point P 4 in figure 1 (b)) will have slight periodic 
changes with time, since the pivoting motion of the wake at the nodal 
line will act to change the perpendicular distance from the point to the 
wake center. The remainder of this section is devoted to mathematical 
statements of the physical picture described in this paragraph. 

Exact model for instantaneous temperature . - The exact solution for 
instantaneous temperature employs a method described in reference 7. The 
molecular wake downstream from a continuously emitting line source is 
considered to be a sheet composed of an infinite number of instantaneous 
line sources. Each of these instantaneous line sources was originally 
introduced at the continuous line source (the heated wire), and the 
degree to which it has diffused into its surroundings is a function only 
of its residence time in the stream, as shown by equation ( 8 ). The 
geometry of the sheet of instantaneous line sources expressed by equation 
(14) is shown diagrammatic ally in figure l(b) and is shown for a partic- 
ular value of the parameter t in the following sketch: 


• ' • 



An important feature shown in the preceding sketch is that the density, 
or proximity, of neighboring instantaneous line sources is decreased when 
the sheet is distorted. This effect is of consequence in later 
discussions. 


The concentration of heat at the fixed point P(X,Y) at time t is 
the summation of contributions of heat from all the instantaneous line 
sources P ! (Ut,y) in the entire sheet. The contribution of heat to 
P(X,Y) by one instantaneous line source at P f (Ut,y) is 


dt 

4 JtCD 


- 55 [ (ut - x)2 + - Y > 2 ] 


ah = 


e 


(15) 
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The total concentration at p(x,Y) is the summation of all contributions: 


Aoo 


h = 


_ sZl 




I e '5 [ (ut - x > 2 + (r - T > 2 ] 


0) 


dt 


(16) 


By employing the functional relations between <x> and y and the resi- 
dence time t (eqs. (5) and (14)) and by changing the independent vari- 
able, equation (16) maybe written as: 


q fl 

4jta 


J_ e '(iS?) ^ {(* A ‘ 4) |s in (*/t ) s in (2 «/t + */t) - -1-1 


y maxJ ( a(*/t) 


(17) 


where y 


max 



(*/) 


An analytic solution to equation (17 ) has not been found; therefore, 
resort must be made to a numerical solution for the exact instantaneous 
temperature. As shown in following sections, a variety of simple alge- 
braic solutions may be used for instantaneous and time-mean temperatures 
under certain sound -fie Id and flow conditions. 

An approximate analytic solution of equation (17 ) has been made for 
points in the mixing region under the following conditions: 

/jj = 4.5, 5.5, 6.5 .. . ^ 

Y = 0 

/t = 0.25, 0.75, 1.25 . . . 


U £ 


4rt2/2cD x 


□. + ^ ) >50 


U2 


) 


(18) 


J 


As shown by the conditions above, the solution is for the maximum instan- 
taneous temperature at the antinodal stations. The solution is derived 
in appendix B, and its final form is 


h » 


q [l 


UK v /4rtoo x 


'1 + 


2vj 

U 2 


where (Ujr is evaluated at the time x/u. 


(for conditions (18) ) 


(19) 
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Approximate model for small deformation of diffu sion wake. - If 
the deformation of the molecular -diffusion -wake center is small, an 
approximate analytic solution can be obtained from the continuous-line- 
source model defined by equation (9). Two conditions are necessary for 
the use of this model: (a) stretching of the molecular wake must be 

small, and (b) the radius of curvature of the deformed molecular -wake 
center must be large with respect to the molecular -wake thickness. The 
following sketch illustrates the necessity for these conditions: 


• • • 

oP 6 . P'(x,y) 






X 


The continuous -line -source model was derived for a perfectly straight 
molecular -wake center plane , so that the density, or spacing, of instan- 
taneous line sources along the wake was constant, or uniform. As shown 
hy the preceding sketch, stretching of the wake causes a departure from 
constant source density. The criterion for condition (a) is that the 


ratio 


Tu 




) 


must be nearly zero, where S is the stretched length 


between nodal points and (u//) is the direct distance between nodal 
points. This ratio is derived in appendix B, and its final form is 



( 20 ) 


where <g is the elliptic integral and is defined as 



( 21 ) 


Equations (20 ) and (21) show that the stretching is a function only of 
the ratio *Jv 2/u. This functional relation is shown in figure 3(a). 
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Consideration of the previous sketch shows that the points Pg and 
Pg receive contributions of different magnitudes from P' and that the 
total contributions received by Pg are greater than those received by 
Pg. Therefore, to approximate the continuous-line-source model, condi- 
tion (b) is necessary. The minimum radius of curvature r c ^ min for any 
time in the cycle is derived in appendix B and has the following value: 


r C,min 


U c 


2 J2 rt/<\/v2 


( 22 ) 


This value occurs when 2rt/x + 2 = rt, 2*, 3*, . . . Condition (b) is 

satisfied when r c min » The minimum radius of curvature is plotted 

in figure 3(b) for a range of values of the parameters #y/v2/u and u//. 
The effect of distance downstream on the ratio r c,min I shown in 

figure 3(c) for a particular value of the coefficient of thermal 
diffusivity a. 


If conditions (a) and (b) are satisfied, then an expression for in- 
stantaneous temperature can be derived from equation (9). Figure 4 is a 
diagram of the molecular -wake position with respect to a fixed point, 
where the distance from the wake center to the fixed point is q. Since 
the radius of curvature is large, the distance T} is approximated by 

T] * (y - Y) sin ^arc cot (23) 


which may he written 



(24) 


Substituting for the lateral distance q in equation (9) and using 
equation (14) for y result in 



( 25 ) 
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which may he reduced to 


2jt 2 / 2 ' 


h * 


q/l 


U/y/4rtO^ 


[”sin sin (zti/'Z + - — — ■ 1 

L \ U / V U / ymaxj 

f 1 + |r sin2 ( 2rt/T + 2 jt/ i)] 


(26) 


If the ratio \/v^/u is very small, then dy/cLx is negligible, 
and the solution becomes 


q/l 2n 2 / 2 cn x 


U^4 jtaj^ 


sm 


(^ t ) sin ( 2ltA + ^ tt ) ' 


7 max 


(27) 


For cases where pc_ is nearly constant over the temperature range 
of interest, equation (2) may be replaced by h = pCp AT, and equations 
(26) and (27) maybe written, respectively. 


-5 [sin (ffi) sin ( 2 «A + ffi) - 


12 


Jmax 


2„ 2 / 2 


“x 


AT e ATq e 


t 


2v^ . 2 

1 + —5— sm' 1 ' 

IT 


^2itA + 2it/^ 


(28) 


AT 


AT 0 e 


2rt 2 / 2 <% 


sin (^j sin ^2 «/t + ^ - 


^max. 


(29) 


where ATq = (q/ 7 )/ (Upcp 4510^.) and ymax = (V2 •/#)/(*/). 

The instantaneous -temperature variations shown in figure 2 were calcu- 
lated from equation (29). 

Another version of the approximate model is described in appendix B. 
In this version, the distance tj shown in figure 4 is determined 
analyt i c ally . 


14 


NACA TN 3760 


Time -Mean Temperature in Mixing Region 


The time-mean temperature at any point in the mixing region may he 
obtained by averaging the instantaneous temperature over one period of 
its cyclic fluctuation: 


AT = / 


' 1 // 


AT dx = 

0 U 


m 


AT d(/x) 


(30) 


where the bar denotes a long-time average. The analysis employs equa- 
tions from the preceding sections , which are written in terms of the 
volume specific enthalpy and can be converted to the temperature differ- 
ence AT = T - T am by means of equation (2). Typical time-mean tempera- 
ture profiles are shown in figure 5 for the case of small deformation of 
the molecular-diffusion wake. An important feature of the time-mean 
temperature profiles is that the mean temperature values are higher near 
the outer edges of the profiles. This is due to the lower transverse 
velocity of the wake center at its extreme excursions. This point is 
illustrated in figure 2 by inspection of the relative magnitudes of areas 
under the temperature variations at the points P-j_, Pg , and P^. Another 


feature of figure 5 is the absence of aerodynamic mixing at the nodal 
points /x/u = 1, 2, . . . 


Calculation of mean temperature for exact model . - Since an analytic 
solution was not found for the instantaneous temperature fluctuations of 
the exact model (eq. (17)), a graphical integration method must be 
employed for evaluation of the mean temperature. Equation (17) may be 
solved graphically for any fixed point P(X,Y) and for a series of 
times x. Then through equation (2), the instantaneous temperature dif- 
ference at the fixed point P(X,Y) may be plotted as a function of the 
time x. The mean temperature difference AT can then be obtained 
graphically. 


Mean temperature for small deformation of diffusion wake . - Calcu- 
lation of time-mean temperature for the case of small wake deformation 
is somewhat simpler than for the exact model, since approximate analytic 
solutions are derived for this case in preceding sections. Instantaneous 
temperature variations, such as those of figure 2, may be plotted from 
equations (28) or (29), and the time-mean temperature may be evaluated 
graphically from such plots. 

Approximate analytic solutions for the mean temperature have been 
obtained from equation (29). The derivations of these solutions and the 
assumptions necessary for their validity are detailed in appendix B. 

The approximate solutions are for the two zones in the mixing region 
shown in the following diagram: 
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where zone I lies within the extreme wake position and zone II lies along 
the extreme wake position. The approximate solutions for zones I and II 
are, respectively, 




q ]l 


2 /\/2p 


c p m l-Ir 


/ 


h/v2 U |sin^n/|j 


1 + 


4F(l-R 2 ) 


jerf [^F (1 + R j] 


erf (1 - R))} - 




1 - R e -F(l+R) 2 . e -F (l-R) 2 
1 + R 


(Zone I) (31) 


where F = 



Close inspection of equation (3l) shows that the time-mean tempera- 
ture AT in zone I is affected by the sound field and stream parameters 

in the following ways when a/f( 1-R) > 2: (a) AT increases as the sound - 

wave frequency rises; (b) AT decreases as the sound intensity rises; and 
(c) AT decreases as the stream velocity rises. Equations (3l) and (32) 
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may be used to construct time-mean temperature profiles such as shorn in 

figure 5 by fairing between zones I and II. As pointed out in the 

RESULTS AND DISCUSSION, the peak value of AT occurs between zones I 
and II, which should be kept in mind when fairing the profile. 


Analytic Procedure for Nonsinusoidal Velocity Fluctuations 

In many sound fields, particularly those of extremely high intensity, 
the velocity fluctuations may depart from sinusoidal form. The analytic 
procedure for this situation is identical to that followed previously in 
the text. If the mathematical expression for the variation of velocity 
fluctuation with time is known, then analytic expressions for the instan- 
taneous and time-mean temperatures might be found. If the velocity 
fluctuation does not follow a known mathematical relation with time, 
then graphical, or numerical, solutions for the instantaneous and time- 
mean temperatures may be made by following the steps in the analytic 
procedure discussed previously. The graphical method is outlined in the 
SUMMARY OF RESULTS AND THEIR PRACTICAL APPLICATION. 


APPARATUS AND EXPERIMENTAL PROCEDURE 
Flow and Sound-Field Measurements 

The apparatus used to produce a standing sound wave transverse to a 
low -turbulence air stream is shown in figure 6. The calming chamber 
minimized the effect of room drafts on the stream velocity and static 
temperature. The stream velocity was controlled by the adjustable -area 
exhaust valve. The standing sound wave was generated by acoustical driver 
units located at opposite ends of the resonant chamber (as shown in fig. 
6). The opposing groups of driver units were connected electrically 
180° out of phase and were powered by a 1400-watt oscillator amplifier. 

The line source of heat was a 0.003-inch-diameter Advance wire continu- 
ously heated by direct current. Calculations of radiation and end con- 
duction showed that these losses were negligible. 

The stream velocity, static temperature, static pressure, transverse 
velocity fluctuations, and sound-wave frequency were measured as described 
in appendix C. These parameters were held constant during the experiment 
at the following values : 
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Stream velocity, 0, ft/sec 17.6 

Ambient static temperature, T^, °F 75 to 82 

Static pressure, in. Hg abs 29.0 to 29.7 

Root -mean-square transverse velocity fluctuation, *Jv 2, ft/sec . . 3.33 

Sound pressure level, db 

Sound-wave frequency, /, cps 104 

Heat addition from line source, q Jl, Btu/ (f t ) (sec ) 0.0308 


Instantaneous and Time -Mean Temperature Measurements 

Two different methods, resistance thermometry and hot-wire anemometry, 
were used to measure the instantaneous temperature downstream from the 
line source of heat. These instantaneous temperature measurements were 
made at several transverse positions in the mixing wake at the axial 
station /x/u = 1.47. 

The sensitive element of the resistance thermometer was a standard 
hot -wire -anemometer probe. The resistance thermometer was essentially 
sensitive only to temperature fluctuations in the airstream. Equations 
describing this instrument’s sensitivity and the compensating circuit 
are given in appendix D. The anemometer probe also was operated as a 
constant-temperature hot-wire anemometer. The hot-wire anemometer was 
sensitive to both temperature and velocity fluctuations in the stream. 
Details of the hot-wire technique and the sensitivity equations are in 
appendix E. 

The mean temperatures at points in the mixing region were measured 
with a thermopile probe as described in appendix F. The mean-temperature 
measurements were made at several transverse positions at the following 
axial stations: /x/u = 0.24, 0.50, 0.75, 1.00, 1.24, 1.47, 1.72, and 
1.96. 


RESULTS AND DISCUSSION 

For the experimental conditions investigated, the stretching of the 
molecular -diffusion wake was small, and the minimum radius of curvature 
was large compared with the wake thickness. Because of apparatus limita- 
tions, experimental verification of the theoretical analysis was possible 
only for the case of small deformation of the diffusion wake. However, 
since large deformations of the diffusion wake are also of practical in- 
terest, this case also is discussed herein on the basis of the theoretical 
analysis alone. 
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Mixing Data and Comparison with Theory for Small 
Deformation of Diffusion Wake 


Flow and sound-field measurements , - Flow and sound-field measure- 
ments were made throughout the region encompassing two nodal stations 
downstream from the line source of heat and transverse distances greater 
than the maximum excursions of the diffusion wake. Measurements with a 
total-static pitot tube and a hot-wire anemometer showed that the axial 
stream velocity was constant throughout the region of interest. The 
hot-wire measurements also showed that the stream turbulence was negli- 
gible and that no vortex street existed downstream from the line source. 


Measurements with the X-wire anemometer showed that the transverse 
velocity fluctuation v in the sound field was nearly sinusoidal and 

that was substantially constant in the region of interest 0 A 

typical oscillogram of the transverse velocity fluctuation is shown in 

figure 7(a), and an axial profile of is shown in figure 7(b). 

The decrease in at the inlet and exit of the resonant chamber was 

due to sound leakage through the openings and sound absorption by the 
chamber walls. 


For the conditions used in the experimental investigation, the 
stretching, radius of curvature, and slope effects had the following 
values : 



0.019 


r C,min/ VW,max - 54.3 


dy 

^ax 


0.19 


These values were obtained from figures 3(a) and (b) for the measured 

values of / \f^ t y U, and o)y ; my (<bx,max is the maximum value of 0 ), 
calculated at the maximum value of X) . The values show that the theo- 
retically predicted wake deformation was negligibly small, as illustrated 
by the following sketch: 


4022 


CL-3 back 


NACA TN 3760 


19 


cvi 

00 

o 



Instantaneous temperatures were calculated from the exact model (eq. (17 )) 
and the approximate model (eq. (29)) for the experimental conditions. 

These calculations showed very close agreement; therefore, the stretching 
and curvature parameters listed above may be considered representative 
for the case of small wake deformation. 

Instantaneous -temperature measurements and comparison with theory . - 
Instantaneous -temperature measurements were made at the second antinodal 
station (/x/u = 1.47) with both the resistance -thermometer and hot-wire- 
anemometer techniques. Typical oscillograms made by these two methods 
are shown in figures 8(a) and (b) , respectively. The measurements were 
made at three transverse positions, Y/y max = 0, 0.49, and 0.98, and show 
good agreement in shape. The hot -wire -anemometer results are discussed 
further in appendix E. 

Instantaneous -temperature variations were calculated from the 
resistance -thermometer oscillograms of figure 8(b). These experimental 
instantaneous -temperature variations are compared in figure 9 with the 
theoretical analysis for small deformations of the diffusion wake as 
expressed by equation (29). Since the transverse velocity fluctuations 
were not precisely sinusoidal, graphical solutions were also made and 
plotted in figure 9. These solutions were based on oscillograms such as 
figure 7(a) and followed the methods outlined by the theoretical analysis. 
The agreement between the experimental data and the theoretical and 
graphical solutions is quite good at the transverse positions Y = 0 and 
Y = 0.17 inch. The discrepancy in figure 9(b) may be due to error in 
the probe position. The peaks in the measured temperature fluctuation 
are closer to /t = 0 than are the theoretical curves, and the measured 
instantaneous temperature is higher than the theoretical value at 
/t = 0. Inspection of figure 2 shows that these two observations are 
consistent with an error in probe position. The asymmetry of the theo- 
retical curves about /t = 0 in figure 9 is due to the fact that /x/u 
was not an exact antinode. Calculation shows that the theoretical 
curve is symmetric about /t = 0.016. 
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Time -mean temperature measurements and comparison with theory . - 
Time-mean temperatures measured with the thermopile probe are compared 
with theoretical analyses in figure 10. The transverse profiles of time- 
mean temperature were measured at a series of axial stations through the 
second nodal station, both with and without the transverse sound field. 
The measured temperatures show good agreement with the equation for 
molecular diffusion (eq c (9)), in the absence of the transverse sound 
field. 


The measured time-mean temperature profiles show good agreement 
with the approximate theoretical solutions, given by equations (3l) and 
(32), except at the axial stations /x/u = 0.75, 1.00, 1.72, and 1.96, 
as shown by figures 10 (c), (d), (g), and (h), respectively. At the 
stations /x/u = 1.00 and 1.96, the measured profiles indicate some 
degree of mixing, contrary to the theoretical prediction of no mixing 
at the nodal stations. The measured profiles at stations /x/u = 0.75 
and 1.72 indicate less mixing than theoretically predicted. These 
discrepancies may be due to probe positioning error or to errors in the 
measurement of stream velocity U and sound-wave frequency / through 
which the theoretical nodal stations /x/u -1,2, . . . are determined. 

The interpolation between the theoretical solutions given by equa- 
tions (31) and (52) is shown by the dotted lines in figure 10. The 
maximum values of AT fall at values of Y/y max less than 1.0. This 

is because at a point Y/y max less than 1.0 the wake passes twice in a 
cycle; and, if the point is near Y/y max = 1.0, the wake passes at a 
velocity less than that at Y/y max = 0. Calculations of time-mean tem- 
perature from the exact model (eqs. (.17) and (30)) showed that the peak 
values of AT fell close to Y/y max = 1.0, as shown in figure 10. 

By closer inspection of figure 10, it can be seen that the agreement 
between theory and experiment is particularly good near the line source 
of heat and at the antinodal stations. It is possible that a random 
angular fluctuation in the stream direction might explain the disagree- 
ment between theory and experiment at far downstream stations where the 
lateral displacements of the wake are small. This possibility is 
strengthened by the observation of low-frequency random fluctuations in 
the thermopile -probe emf readings throughout the course of the experiment. 

Graphical solutions, based on the velocity fluctuation oscillograms 
(fig. 7(a)), were made for the axial stations /x/u = 0.75 and 1.47 as 
shown in figures 10(c) and (f), respectively. These graphical solutions 
are in fairly good agreement with the analytic solutions based on a 
sinusoidal velocity wave form, so the nonsinusoidal wave form was not 
the reason for disagreement between theory and experiment. 
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The discrepancies discussed appear to he caused hy imperfect experi- 
mental conditions and techniques. For this reason, the experimental data 
are considered to be satisfactory substantiation of the theoretical 
analysis for small deformation of the diffusion wake. 


Theoretical Effect of Large Diffusion-Wake Deformations on Mixing 

In order to illustrate the effect of wake deformation on the instan- 
taneous and time-mean aerodynamic mixing, a comparison was made between 
the exact model and solutions applying in the case of negligible wake 
deformation. The exact model (eq. (17)) accounts fully for stretching 
and curvature of the wake. The solutions for negligible wake deformation 
(eqs. (29), (28), and (B40) ) are predicated on the assumptions of negli- 
gible wake stretching and negligible diffusion in the stream (axial) 
direction. In addition, equation (29) assumes that the wake center is 
nearly parallel to the stream direction (i.e., small slope), and equation 
(28) assumes that the wake center has negligible curvature. Equation 
(B40) accounts for the geometric curvature of the wake and is derived in 
appendix B. 

In order to determine the effect of appreciable wake deformation on 
the mixing parameters AT and AT, solutions from the exact and approxi- 
mate models were calculated for the following conditions : 

= 16.7 fps (161 db) 

/ as 520 cps 
U = 17.6 fps 

q/l a* 0. 0308 Btu/ (ft ) (sec ) 

/X/U = 7.5 

These conditions were chosen so that the seventh antinode (/x/u =7.5) 
fell at the same axial station as the second antinode in the experiment. 
This kept the width of the molecular wake the same as in the experiment. 
The conditions above are also such that y max is the same as in the 

experiment so that a graph of ZiT against A would show local mixing 
in space. The stretching and curvature parameters for the conditions 
above are 


- 0.36 

r n . / -\A> = 2. 

Cornin' v max 


22 


NACA TN 3760 


Effect on instantaneous temperature . - Instantaneous -temperature 
differences AT were calculated from the exact model (eq. (17)), for 
the conditions above , at a series of transverse positions. These instan- 
taneous temperatures are shown by the solid lines in figures 11 (a) to 
(d). For comparison with the exact model, instantaneous -temper ature 
differences were calculated for the same conditions from the simplest 
solution of the small-deformation model as given by equation (29). The 
small-deformation-model solution is shown by the short-dash lines in 
figures 11(a) to (d). This comparison of exact and approximate models 
clearly shows that a substantial degree of instantaneous, or local, 
mixing occurs when the wake deformations are large. For all the trans- 
verse positions shown, the peak temperature of the instantaneous time 
profile is substantially lower when the wake is deformed. This decrease 
in peak temperature is due to the stretching of the wake and concurrent 
lower density of instantaneous sources, as discussed in the THEORETICAL 
ANALYSIS. The decrease is largest for the transverse position Y = 0, 
since the source density is lowest for the shape assumed by the wake as 
it passes through Y = 0 at an antinodal station. This point can be 
verified from inspection of the sketches in the THEORETICAL ANALYSIS 
section and of figure 4. 

In order to illustrate the effect of curvature, equation (28) is 
plotted in figure 11 as a broken line. Equation (28) accounts for the 
slope of the wake center, thereby being a more accurate representation 
of the small- deformat ion model. In addition, a precise representation 
of the small-deformation model is shown in figure 11 by the long-dash 
lines. This solution (eq. (B40)) accounts fully for the geometric curva- 
ture of the wake center, as shown in figure 4, and is derived in appendix 
B. The effect of curvature on the exact model is illustrated in figure 
11(c). By inspection of figure 1(b), it is seen that a fixed point at 
y/y max =0 in the antinodal plane receives significant contributions 

from a longer arc length of the wake when the wake is at f x = 0 than 
when the wake is at fx = 0.15. Because of this curvature effect on 
contributions to Y/y max = 0.8, the exact model shows a correct tempera- 
ture at fx = 0 which is higher than any of the solutions to the small- 
deformation model. This behavior of the exact model in comparison to 
that of the small-deformation model is quite reasonable in view of the 
fact that the small-deformation model implicitly assumes a flat, uncurved, 
center shape. 

The peak-instantaneous -temperature difference calculated from 
equation (19) is also shown in figure ll(a). The comparison of equation 
(19) with equation (17) is quite good, indicating that equation (19) may 
be used for approximate calculation of the peak instantaneous temperature 
for the case of large wake deformation. Algebraic comparison of equa- 
tions (19) and (27) shows that the peak temperature at the antinodal 
stations and Y = 0 is diminished approximately by the factor 



when wake deformation becomes appreciable 
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In summary, large wake deformations substantially lower the peak 
values of the instantaneous temperature and appreciably change the shape 
of the temperature-time profile. 

Effect on time -mean temperature . - The effect of wake stretching 
and curvature on the time-mean temperature is shown in figure 12, where 
theoretical temperature differences Zfl? are plotted against the trans- 
verse position Y/y max . The values of mean temperature difference were 

calculated from equation (30 ) for each of the four solutions shown in 
figures 11(a) to (d). 

As discussed in the preceding section, equation (B40) is the precise 
representation of temperature fluctuations for the case of negligible 
wake deformation. The mean temperatures calculated from the exact model 
(eqs. (17 ) and (30)) are considerably lower than those calculated from 
the negligible deformation model (eqs. (B40) and (30)). This comparison 
clearly shows the effect of stretching and curvature on the time -mean 
mixing . 

The instantaneous -temperature variations calculated from the negli- 
gible deformation model (eq. (29)) are very different from those calcu- 
lated from the exact model (eq. (17)), as shown in figure 11. For this 
reason, the close agreement between the mean temperature differences 
calculated from these two models and shown in figure 12 is fortuitous. 

In summary, large wake deformations substantially lower the time- 
mean temperatures throughout the transverse profiles in the mixing region. 


SUMMARY OF RESULTS AND THEIR PRACTICAL APPLICATION 


In order to facilitate the engineering application of the informa- 
tion presented in this report, the results are summarized in this section. 
The information is divided into four practical problem cases : for sound 

fields with a sinusoidal wave form (l) for small wake deformation and 
(2) for large wake deformation, and for sound fields with a nonsinusoidal 
wave form (3) for small wake deformation and (4) for large wake deforma- 
tion. The criteria for small wake deformation are low values of the 


stretching parameter — — ^ can evaluated from figure 3(a), 


C,min 


/V a \n; 


ax* 


and high values of the curvature parameter 
radius of curvature ^ ^ 

spreading coefficient o> can be evaluated from equation (5). 


where the 


r P can be evaluated from figure 3(b) and the 

o ,min 

if 


1 C ,min / V^max 
negligible. 


is greater than 50, then the effect of curvature is 
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Another criterion for small wake deformation is the comparison of 
peak instantaneous temperatures at the antinodal stations and Y = 0 
calculated from equations (19) and (27), (Note conditions imposed on 
eq. (19) in the THEORETICAL ANALYSIS section.) This comparison shows 
that if (2v 2 /u 2 ) ~ 0, then the peak temperatures calculated from these 
two equations are nearly equal and the wake deformation is small. 

In the sections following, equations and procedures are listed for 
each of the four cases discussed previously. 


Sound Fields with Sinusoidal Wave Form 

Small wake deformation . - The conditions necessary for small wake 
deformation are 



r C,min > > 


^%ax 


A /^2 

U 


* 0 


Recommended equations: 


ZXT(X,Y,t) 




fl/Y 

V2 V y2_ 


2 


(29) 
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Large wake deformation: 



' ^7t|( n/t " n/ |) + [sin(jc/t)sin(2itA + */t) 

1 l L -\/2 yvZ. 


( 17 ) 


From equations (17) and (2), plot AT(X,Y,iO against T. Then: 

^l// 

AT(X,Y) = / AT dT 


(30) 


Sound Fields with Nonsinusoidal Wave Form 
Small wake deformation . - Conditions necessary: 

0 


Recommended procedure: 
From: 


S - (U//) 

m) 


C,min >> V°^iax 

V7 „ _ 


u 
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perform graphical integration to obtain: 



Using equation (12), plot y(x/t): 
Fraction of 



For a particular X and Y, by using equations (9) and (2) and method 
shown in figure 4 and described by equation (B40), find: 
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Large wake deformation, - Recommended procedure: 


From: 



perform graphical integration to obtain: 



Using equation (12), plot y (x,r ) 


Fraction of 



*(X,Y/0 = I £ 


1 


- TT Q Ut - X ) 2 + b - Y ) 2 ] 

dt 


(16) 
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For a particular X and Y, from equations (16) and (2), plot: 



/ 


For a particular X, by graphical integration, find: 



CONCLUSIONS 

From an investigation of the aerodynamic mixing by a standing sound 
wave downstream from a continuous line source of heat, the following con- 
clusions are drawn: 

1. The velocity fluctuations in a periodic sound field contribute 
to aerodynamic mixing by periodically displacing the diffusion wake in a 
manner similar to the displacements of a flag waving in a harmonic mode. 
Important characterizing parameters of the wake displacement in a 
sinusoidal-wave -form sound field are: 

a. Nodal points of the wake displacement exist at the 
stations /x/u = 1 , 2, . . . 

b. The maximum wake displacement occurs at the antinodal 

stations and is equal to (V2 x rms transverse velocity fluctuation) 
divided by (n x sound-wave frequency). 

2. When the stretching of the diffusion wake is small and its 
radius of curvature is much greater than its thickness, then the wake 
deformation is small, and the sound field has the following mixing 
effects : 
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a. In the time-mean sense, concentrations are greatly- 
reduced, except at the nodal stations. For example, near 

the line source the time-mean temperature difference is reduced 
by more than an order of magnitude in the presence of a 147- 
decibel, 104-cycle-per-second sound field. 

b. In the instantaneous , or local, sense the sound field 
does not contribute to mixing. 

3. When the sound field causes appreciable deformation of the dif- 
fusion wake by stretching and bending, then the sound field contributes 
to mixing in both the time-mean and instantaneous senses. For example, 
the peak instantaneous temperature difference in a 161-decibel, 520-cycle- 
per-second sound field was 61 percent of the peak instantaneous - 
temperature difference in a 147-decibel, 104-cycle-per-second sound 
field; and the mean temperature was as low as 66 percent of that in the 
147-decibel field. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, May 9, 1956 


APPENDIX A 


SYMBOLS 


term in King's eq. 

constant in King's eq. 

term in King's eq. 

constant in King's eq. 

capacitance 

concentration of mass 

specific heat at constant pressure 

coefficient of molecular diffusion of mass 

hot-wire diameter 

transverse displacement of fluid line elements 

sound pressure level in decibels = 20 l°8]_o 
jp/ (0.0002 dynes/sq cm)) 

direct -current voltage 

complete elliptic integral of second kind 
alternating-current voltage 


error function 



factor equal to 



constant 


frequency 


electronic gain 
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g 

E/l 

h 

I 

k 

l 

M 

m 

P 

P* 


P 

qjl 

R 


r C 

S 

T 

AT 

AT 

t 

U 

u 

v 


constant 

strength of instantaneous heat line source per unit length 
volume specific enthalpy (enthalpy per unit volume) 
direct current 
thermal conductivity 
length 

time constant 

constant in expression for film temperature 
fixed point in mixing region 

position of instantaneous line sources in molecular-diffusion 
wake 

sound pressure fluctuation 

rate of heat release from line source per unit length 

factor equal to (Yjt /) 
radius of curvature of deformed molecular-diffusion-wake center 
stretched internodal length of molecular-diffusion-wake center 
temperature, °R 

temperature difference, T - T^, °F 

time-mean temperature difference 
residence time 
stream velocity 
axial velocity fluctuation 

transverse velocity fluctuation, (i.e., particle velocity in 
sound wave) 
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X,Y 

position of probe, or fixed point 

x,y,z 

Cartesian coordinates fixed in space 

y 

position of molecular-diffusion-wake center with respect to 
flow axis plane 

Ymax 

maximum position of molecular -diffusion -wake center with 
respect to flow axis plane = yv 2 )/(jt/) 

a 

coefficient of thermal diffusivity, k/pc p 

P 

variables of integration 

r 

Gamma function 

6 

temperature fluctuation 


Cartesian coordinates moving with the stream 

e 

an angle 

p 

density 

a 

thermal coefficient of resistance 

T 

time associated with the standing sound wave 

♦ 

angle of X-wire array to stream direction 


electric resistance 

(D 

molecular spreading coefficient 


molecular spreading coefficient evaluated at x/U 


Subscripts: 


A 

amplifier input 

am 

ambient 

c 

compensator 

f 

gas film 
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i summation convention 

L electric cables 

max maximum 

min minimum 

o output 

w referring to hot wire 

0 reference 

1 any particular value for the variable 

Superscript : 


mean 
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APPENDIX B 

DERIVATION OF EQUATIONS APPEARING IN THEORETICAL ANALYSIS 
Approximate Analytic Solution of Exact 
Instantaneous -Temperature Model 

Equation (17) can be solved analytically for the special conditions: 

ic/t - jc/ | <0.02 */ | (Bl) 


/ = 4.5, 5.5, 6.5 . . . 

(B2) 


Y = 0 

(B3) 


/ t = 0.25, 0.75, 1.25 . . . 

(B4) 

- 

U 2 L + >50 

(B5) 

- 

4jt 2 / 2 co x \ U 2 / 



Letting P = n/t - rt/ ^ and using equation (13) for y in equation (16), 
result in 



It should be noted here that the application of condition (Bl) limits the 
range of integration to the major values of the integrand. With condi- 
tions (Bl), (B3), and (B4), an approximate form is 


h 


q/7 

4* 2 /a> x 


r >ao 



u 


2 


4jt 2 / 2 a>x 




(B7) 
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By trignometric identity and conditions {B2) and (B4), 


oo 

CV3 

O 


h * 


q Jl 


4rt 2 /u^ 




4jt 2 / 2 CD x 


3 2 + (i sin 23 Y 

ATT^ ' / 


4U 


d3 


/ x 
-V u 


If the angle p is limited to the conditions 

|B| <0.25 

then the sine of 23 * can be replaced by 23: 


(B8) 


(B9) 


r'ao 


qA 


4rt 2 /' 


oo.. 


4^ + U 2 ' ^ 


d3 




(B10) 


The imposition of condition (B9) is justified by condition (B5), which 
drives the integrand to small values at the limits of condition (B9). 
By change of variable, equation (B10) may be written 



By conditions (B2) and (B5), and by algebraic reduction, 

h « — — — — (for conditions (Bl) through (B5)) (B12) 

U V^ y x + 2v£ 
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Calculation of Stretching of Molecular-Diffusion Wake 


The stretched length S of the molecular -diffusion wake between 
adjacent nodal points in the mixing region can be calculated from the 
equation 

•■fRC 

The derivative dy/dx can be determined from equation (12) in the THEO- 
RETICAL ANALYSIS section, so that the stretched length is 


1/2 


dx 


(B13) 
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Since sin 2 0 = cos 2 (p - and cos 2 0 = 1 - sin 2 0, equation (B14) may be 
written 



By the change of independent variable cp = 2 jc/t + 2ji / - p and by noting 

the periodicity of the integrand, equation (B15) becomes 


S 



vJO 


1/2 



dtp (B16) 


The integral is the 


complete elliptic 


integral of the second kind; 



(B17) 
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so that the stretched length between adjacent nodal points is 


S 



+ 



g 


From equation (B18) a stretching parameter may be written 
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(B18) 
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Calculation of Radius of Curvature of Deformed 
Molecular -Diffusion Wake 
The formula for radius of curvature is 



(B20) 


Taking the first and second derivatives of equation (ll) and substituting 
into equation (B20) result in 


[«.(** . wj)| 


where the vertical bars indicate absolute values. The minimum radius of 
curvature r c ^ min occurs when 2rt/c + 2«/ ^ . . . , and its 

valu^ is 


r C,min 
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Approximate Solutions of Mean-Temperature Integral 
for Negligible Wake Deformation 

As shown in t he THEORETICAL ANALYSIS section, an expression for the 
mean temperature AT at any point in the mixing region is the integral 
equation formed from equations (29) and (30): 



where 

w - (v* ipy™ 

Although an analytic solution could not he found for equation (B23), ap- 
proximate solutions were found for the two regions shown in the following 
diagram: 



The conditions for approximate solutions in these two zones are: 
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Approximate solution for zone I. - By a change of variable, equation 
(B23) may be transformed to 


— at 0 


p2rt + */ £ 


AT = 


2jt 


U 


a -F(sin x - R)‘ 


dx 




v 2 sin 1 


(B24) 




/• r X * ^ \ U / ^ Yn/ 

where X = 2jt/x + Jt/ — , F = 9 > a^d R = ~ r==. /rr/xN ’ 

u 2rt 2 / 2 (J0 Y2V^ sin ( u) 

Since the integral limits represent one cycle, equation (B24) may be 
written as 


_ at 0 

AT = — - 


‘jt/ 2 
-it/ 2 


-F(sin x - r) 2 dK 


By another change in variable, equation (B25) becomes 


(B25) 


— t — AT 0 
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AT = 


jt 


-F(e - R)‘ 


de 


(B26) 


_! Vl - 

where e = sin x. With a third change in variable, equation (B26) becomes 

' - 1-R 

^-Fv 2 

dv (B27 ) 


— AT 0 
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* 1 - r2 ) - ( y2 + 2Rv ) 
1 -1-R 


where v = e - R. Expanding the integrand denominator in a binomial 
series in (l - R 2 ) and (v2 + 2Rv), evaluating at the integral limits, 
and neglecting higher order terms give 


F^ 1 + + erf[VF(l-R)]| - 
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Noting that 


^0 = 


q fl 


pCpU 


equation (B28) may be further reduced to 



Approximate solution for zone II . - When Y 
then R = -1, so that equation (B27 ) becomes 



AT « 




(B30) 


Expanding the integrand denominator in a binomial series and neglecting 
higher order terms , 


AT * 





e~ 


dv 


By a change in variable, equation (B3l) becomes 


(B31) 


AT 



(B32) 


As shown in reference 8, when ^ 2 . fV 4 >1, equation (B32) is approximated 
by 
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Solution for Instantaneous Temperature for Small Wake 
Deformation Allowing for Wake Curvature 

The THEORETICAL ANALYSIS section presents two approximate solutions 
for small deformations of the molecular-diffusion wake. Figure 4 shows 
the geometric assumptions used in deriving these equations. In equation 
(26), the slope of the wake center plane was assumed constant over the 
small range of tj which is of interest. Derivation of equation (27) as- 
sumed further that the slope was essentially zero, thereby making it 
possible to neglect the contribution of the axial displacement to t). 

The following derivation considers the curvature of the wake center 
plane. The position of the wake at any particular fraction of a cycle 
A i is given by 


y = sin (it) sin ( 2 * 710 . + (B34) 

There is a point (x-^,y^) satisfying equation (B34) which is a minimal 

distance T] from the probe at (X,Y). This minimal distance condition 
requires that t] be measured normal to the center plane curve, as shown 
in figure 4. Therefore, the equation for the slope and end points of the 
line segment t] must be given by 



(B35) 


where 





Combination of (B35) and (B36) gives 


(B36) 


- Y = 


U 


(X - x x ) 


s in ( 2 * /x-^ + — n~) 


(B37) 


Equations (B34) and (B37) define all that is needed to calculate the 
minimal distance rj for any time T . Direct substitution of q into 
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equations (9) and (2) gives the instantaneous temperature at any position 
within the mixing region. The actual procedure used to obtain figure 11 
is outlined in the following paragraphs. 

The solution of equations (B34) and (B37 ) simultaneously for Xq is 


sin 





(X - Xq) 



2rt/xq 


«/ 




U 


= 0 


(B38) 


For assigned values of Xq near the probe-position plane X, equation 
(B38) can be solved for Aq by successive trials. The corresponding 
value of yq is obtained from equation (B34), and the calculation of q 
at this particular ] iq can be made using the following equation. 

T|q = [(X - Xq ) 2 + (Y - yq) 2 ] l/2 ( B39 ) 

The subscript 1 indicates that the distance corresponds to a particular 

Aq. 

The instantaneous temperature AT^_ a.t this fraction of a cycle may 
be made by using the following simplification of equations (9) and (2): 

ATj. - «T 0 CB40) 


where AT 0 = (q/l)/(V^ P c p U )> and ^1 is S iven by equation (B39) . 

Equation (B40) is restricted by the two conditions discussed in "Approx- 
imate model for small deformation of diffusion wake." If these conditions 
are satisfied, this analysis is a better approximation of the exact solu- 
tion (eq. (17)) than either equations (28) or (29). However, it is ob- 
viously a far more tedious procedure. 
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APPENDIX C 


FLOW AND SOUND -FIELD MEASUREMENTS 
Stream Velocity 

The stream velocity was measured ty the vortex-shedding technique 
described in reference 9. The frequency of vortex shedding from a 0.250- 
inch-diameter cylinder was determined by forming a Lissajou's figure 
from a hot-wire anemometer signal. The stream velocity was then calcu- 
lated from the frequency of shedding by means of the Strouhal number. 


Stream Static Temperature and Pressure 

The stream temperature was measured with a mercury thermometer at 
the calming-chamber inlet. The stream pressure was measured with a ba- 
rometer and differential water manometer. 


Transverse Velocity Fluctuation 

The transverse velocity fluctuation v was measured with an X-wire 
anemometer. Since the velocity antinode and pressure-node planes of the 
standing sound wave passed through the line source and stream centerline, 
the velocity fluctuations were at a maximum and the pressure fluctuations 
were negligible in the mixing region. Since both turbulent velocity 
fluctuations and pressure fluctuations were negligible, the X-wire ane- 
mometer could be used to measure the transverse velocity fluctuation. 

The operation and design of the X-wire anemometer equipment has been 
described in references 1 and 10, and will be only outlined here. The 
X-wire probe and circuit diagram were identical to those described in 
reference 10. The sensitivity of the X-wire is given by the equation 

v = gIZI/a^T e (cl) 

where e is the alternating -volt age difference between the two hot wires, 
G is the amplifier gain, and AE/A\|r is the slope of the X-wire calibra- 
tion curve. A typical plot of the direct -current -voltage difference be- 
tween the hot wires E and the relative angle between the X-array and 
the stream direction i|r is shown in reference 1. 


The root -mean-square value of the transverse velocity fluctuation 

was measured with a true root -mean-square voltmeter. The sound- 
wave velocity v(t) was also recorded in the form of oscillograms. 
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APPENDIX D 


INSTANTANEOUS -TEMPERATURE MEASUREMENTS WITH RESISTANCE THERMOMETER 

A standard single-wire anemometer probe was used as the sensitive 
element of the resistance-thermometer equipment and is shown in figure 
13(a). The resistance thermometer circuit design was based on unpublished 
instrument research by Mr. E. Carlson of the Lewis laboratory. The circuit 
diagram for the resistance thermometer is shown in figure 13(b). With 
the Tf run-calibrate switch" in the calibrate position, a step-function 
current could be impressed on the resistance thermometer. A typical time 
response of the resistance element to this step-current is shown in fig- 
ure 14(a). Since the resistance element was known to have an exponential 
response curve, its time constant could be measured directly from oscil- 
lograms as shown in figure 14(b). 

With the "run-calibrate switch" in the run position, the circuit in 
figure 13(b) could be used as a compensated resistance thermometer. The 
sensitivity of the instrument to temperature fluctuations is given by 
the expression 


'•'-ami 

[0.1404 (2 am + 2 l ) + 1] 

1 e 

31.62 

^c^-w^am j 

f ^am 1 

L 1 + CT ( T am " 

T o)J 


(EL) 


where 6 is the temperature fluctuation, T flm is the ambient air tem- 
perature, 2^ is the wire resistance at the temperature T^, 2^ is 
the probe lead resistance, G is the gain of the amplifier, I w is the 

wire current, a is the thermal coefficient of resistance of the wire 
material, e is the amplifier output voltage, and 31.62 G c is the gain 

of the resistance thermometer circuit where G c is given by the 
expression 


G c = 



(D2) 


In the present experiment, the circuit components were adjusted so that 
the term (2it/2 c C)^ was small compared to the first term in the denomi- 
nator. The gain G c could then be approximated by the expression 


I 


I 
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In operation, the time constant M was determined from an oscillo- 
gram of the time response of the resistance element to the step current. 
The resistances 2 C and Qq were then adjusted so that 

2 c C = M (D4) 

and 

^ = 14 ( D5 ) 

“0 

With the condition (D5), the over-all circuit frequency response was 
flat within 3 percent from 20 to 500 cycles per second. With the cor- 
rect compensation, the circuit was then used to obtain oscillograms of 
the instantaneous -temperature fluctuation at a number of points in the 
mixing region. 
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APPENDIX E 


INSTANTANEOUS -TEMPERATURE MEASUREMENTS WITH HOT-WIRE ANEMOMETER 

As pointed out in reference 11, the hot-wire anemometer theoreti- 
caipy can he used to measure temperature fluctuations simultaneously 
with velocity fluctuations. In the present experiment, the axial veloc- 
ity fluctuations were negligible, so that a single hot wire placed normal 
to the flow direction should be sensitive only to temperature fluctuations. 
However, for clarity and future reference, the hot-wire sensitivity equa- 
tion is derived herein for both temperature and velocity fluctuations. 

The derivation is made for constant -temperature operation of the hot wire. 


Derivation of Sensitivity Equation 


In reference 11, King’s equation is written: 



I 2 2 = (a + b Vu)(2 - SqJ 

(El) 

where I and 
respectively. 

are functions 

2 are the hot-wire current and operating 
The factors A and B and the cold-wire 

of temperature as follows: 

resistance, 
resistance 2 am 


A = #- k 

V 

(E2) 


B - v 

(E3) 


“am " “ot 1 + °( T am - V] 

(E4) 


where a and b are constants, l is the wire length, is the wire 
diameter, a is the thermal coefficient of resistance of the wire mate- 
rial, is the wire resistance at the reference temperature Tq, ^am 

is the wire resistance at the ambient temperature T am , k is the fluid 
thermal conductivity, c p is the fluid isobaric specific heat, and p is 
the fluid density. 

In the temperature range from 500° to 1000° R, the thermal conduc- 
tivity of air is closely approximated by 

am 


k = fT 


0.851 


(E5) 
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where f is a constant , and the exponent 0.851 was evaluated from the 
data of reference 12. Since the specific heat Cp has negligible change 

in the temperature range 500° to 1000° equations (E2) and (E3) can be 
written as 

0,851 <e6) 

B = (^V*^e)[ T - m(T - t ^O ' 0,0745 (E7) 

where jjF - m(T - T njn >] is the film temperature appropriate to the heat 

transfer from the hot wire at temperature T to the fluid stream at 
temperature T am and m is an unknown constant . 

With A, 3, and 2^ given in terms of the temperature equa- 

tion (El) can be differentiated to obtain the combined temperature and 
velocity sensitivity of a constant-temperature hot-wire anemometer (2, 
and therefore T, held constant): 


212 dl = J- [( A + B VU)(8 - “am>] ^am + |j [(* + B Vu)(“ - “am)] 
u± am 

(E8) 

By noting that the voltage fluctuation across the hot wire is e = 2dl 
and defining the temperature and velocity fluctuations as & = dT^ and 

u = dU, respectively, equation (E8) becomes 


2Ie = 


(A + BVu) 5^— (2 - 2^) + 

u± am 


(2 ' 2ajn) + ^ 5t^)_ 


8 + 


> - 2 am) E l" 

2 Vu - 


(E9) 


The various partial derivatives are 


d 
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a£, 


am 


(2 - 2 am ) - -°2 0 - - 1 + a(T am - T 0 ) 


(E10) 


48 


NACA TN 3760 


$>A 



(0.851 



- “(T - Tgj] 


-0.149 


0.851 mA 
[t - m(T - T m )] 

(Ell) 


dB _ 0.0745 mB 

^am " ' [T - m(T - Tj] 


(E12) 


By assuming negligible natural- convection effects, the mean values A 
and B can be written as 


A = 


lo 2 


2 - 

am 


B = 


(I 2 - I§) 

VU(2- 8 m ) 


(E13) 


(E14) 


where Iq is the wire current at the stream temperature and pressure, 
but with zero velocity. 

Substitution of equations (E10) through (E14) into equation (E9) re- 
sults in the desired sensitivity equation 



_ am 


(E15) 


Hot -Wire -Anemometer Circuitry and Operation 

The hot -wire-anemometer probe used for the temperature fluctuation 
measurements is shown in figure 15(a), and the circuit diagram was the 
same as that used in reference 10. Since the velocity fluctuations in 
the axial direction had negligible effect on a single hot wire positioned 
normal to the stream direction, equation (E15) could be modified to re- 
late the oscilloscope display voltage to the stream temperature fluctua- 
tions as follows: 
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where 2G is the gain of the bridge-amplifier combination, and e 0 is 
the voltage displayed on the oscilloscope. 


Hot -Wire -Anemometer Measurements 

The sinusoidal wave to be seen at the base of the hot-wire-anemometer 
oscillograms in figure 7 (a) is due to the influence of the transverse ve- 
locity fluctuation v on the hot-wire heat loss. As shown above, the 
sensitivity of the hot-wire anemometer to temperature fluctuations depends 
on the film temperature Tp = T w - m(T w - T^) . When the constant m is 

set equal to 1.0 (so that Tp = T^), and when the oscillograms are cor- 
rected for the effect of v, then the peak temperature difference ATq 

measured by hot-wire anemometry shows good agreement with the peak temper- 
ature AT n predicted by equation (7). The equivalence of Tf and T _ 

^ i am 

is in agreement with previous results reported in reference 13, where the 
term A (eq. (E2)) was found to he a function of the ambient fluid 
temperature. 

In summary, the hot-wire-anemometer measurements substantiate those 
of the resistance thermometer when the film temperature was taken equal 
to the ambient temperature. Further research is necessary before the 
equivalence of film and ambient temperature is fully established. 
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APPENDIX F 


MEAN -TEMPERATURE MEASUREMENTS 

The mean-temperature measurements were made with an iron-constantan 
thermopile probe shown in figure 15(a). The thermo junctions were care- 
fully alined to fall on a line parallel to the line source of heat. The 
circuit diagram is shown in figure 15(b). The time constant of the ther- 
mocouples was high enough to ensure true averaging of temperature fluc- 
tuations, so that the emf read with the potentiometer was in direct 
proportion to the mean temperature difference AT. 
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(a) Typical temperature profiles in molecular- diffusion wake at an instant. 
Figure 1. - Position of molecular-diffusion wake in mixing region. 
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(b) Position of molecular-diffusion-wake center at various times through one cycle 
of transverse -velocity fluctuation. 

Figure 1. - Concluded. Position of molecular -diffusion wake in mixing region. 


NACA TN 3760 


Dimensionless temperature difference, AT/AT q 
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Fractions of transverse-velocity- 
fluctuation cycle, / t 

(c) At point P 3 in figure l(h) . 

Figure 2. - Typical theoretical instantaneous - 
temperature - time variations in mixing region. 
Dimensionless axial position, 1.5; coefficient of 
thermal diffusivity, 2.37xl0 -4 square feet per sec- 
ond; sound-wave frequency, 104 cycles per second; 
root -mean- square transverse velocity fluctuation, 
3.33 feet per second; rate of heat release from line 
source per unit length, 0.0308 Btu per second per 
foot. 
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Ratio of rms transverse velocity to 

a/v2 

stream velocity , v 


(b) Minimum radius of curvature of molecular - 
diffusion wake. 


Figure 3. - Continued. Factors characterizing 
deformation of molecular -diffusion wake. 
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(c) Curvature parameter; coefficient of thermal diffusivity, 
2.37X10” square feet per second. 


Figure 3. - Concluded. Factors characterizing deformation of 
molecular-diffusion wake. 
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Figure 4. - Local diagram of position of molecular-diffusion-wake 
center at various times. 
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Figure 5. - Typical time-mean temperature profiles in mixing region for small 
deformation of the molecular- diffusion wake. 
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Figure 6. - Flow and sound-field apparatus. 
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Proportional to time t 


(a) Typical oscillogram of transverse velocity fluctuation. 



(b) Axial profile of velocity fluctuation in transverse standing 
sound wave. Reference root -mean- square transverse velocity 
fluctuation at x = 2 inches. 

Figure 7. - Transverse -velocity-f luctuation data. 
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(a) Hot-wire-anemometer data. 
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y max 
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Proportional to time T 

(b) Resistance- thermometer data. 

Figure 8. - Instantaneous-temperature-difference oscillograms, 
inches or /X = 1.47. 
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-.5 — • 4 -.3 -.2 -.1 0 .1 .2 .3 .4 .5 

Fraction of transverse-velocity-fluctuation cycle, /t 

(a) Lateral position, Y = 0 inch. 

Figure 9. - Instantaneous temperature-time variation at various lateral positions. 
X = 2.98 inches or /X/U = 1.47. 
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-.5 -.4 -.3 -.2 -.1 0 .1 .2 .3 .4 .5 

Fraction of transverse-velocity-fluctuation cycle, /t 

(b) Lateral position, Y = 0.085 inch or Y/y max = 0.49. 

Figure 9. - Continued. Instantaneous temperature-time variation at various lateral 
positions. X = 2.98 inches or /X/U =* 1.47. 
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-.5 -.4 -.3 -.2 -.1 0 .1 .2 .3 .4 .5 

Fraction of transverse-velocity fluctuation cycle, /t 

(c) Lateral position, Y =0.17 inch or Y/y max = 0.98. 

Figure 9. - Concluded. Instantaneous temperature-time variation at various lateral 
positions. X = 2.98 inches or / X/U = 1.47. 
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(a) Axial station, X = 0.49 inch or = 0.24; (b) Axial station, X = 1.02 inches or ^ = 0.50; 

y m ax “ °- 119 inch ‘ U y m ax “0.173 inch. 

Figure 10. - Comparison of theoretical and experimental time-mean temperature- difference 
profiles at various axial stations. 
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(c) Axial station, X = 1.52 Inches or — - 0.75; 
y mav * 0.123 Inch. 


(d) Axial station, X = 2.03 Inches or 

^ = 1 . 00 . 


Figure 10. - Continued. Comparison of theoretical and experimental time-mean temperature-difference 
profiles at various axial stations. 
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(e) Axial station, X = 2.51 inches or = 1.24; 
y = 0.119. u 


/ y 

(f) Axial station, X = 2.98 inches or = 1.47; 

y =0.172 inch, 
max 


Figure 10. - Continued. Comparison of theoretical and experimental time-mean temperature- 
difference profiles at various axial stations. 
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(g) Axial station, X « 3.49 inches or = 1.72; (h) Axial station, X « 3.98 inches or = 1.96. 

y max " °* 122 lnch ‘ 


Figure 10. - Concluded. Comparison of theoretical and experimental time-mean temperature- 
difference profiles at various axial stations. 
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(b) Temperature- time variation at Y/y max =« 0.4. 

Figure 11. - Comparison of exact and approximate solutions for instantaneous tenmerature. Dimension- 
less axial position, 7.5; coefficient of thermal diffusivity, 2.37X10 -4 square feet per second; 
stream-velocity - frequency ratio, 0.0338 foot; ratio of root-mean-square velocity fluctuation to 
stream velocity, 0.947; rate of heat release from line source per unit length, 0.0308 Btu per second 
per foot. 
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(c) Temperature- time variation at Y/y ma x " °* 8 * 



Figure 11. - Concluded. Comparison of exact and approximate solutions for instantaneous temperature. 
Dimensionless axial position, 7.5; coefficient of thermal diffusivity, 2.37X10" 4 square feet per 
second; stream- velocity - frequency ratio, 0.0338 foot; ratio of root-mean-square velocity fluctua- 
tion to stream velocity, 0.947; rate of heat release from line source per unit length, 0.0308 Btu 
per second per foot. 
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Figure 12. - Comparison of exact and approximate solutions for mean temperature. Dimensionless axial position, 7.5; 
coefficient of thermal diffusivity , 2.37X10" 4 square feet per second; stream-velocity - frequency ratio, 0.0338 
foot; ratio of root -me an- square velocity fluctuation to stream velocity, 0.947; rate of heat release from line 
source per unit length, 0.0308 Btu'per second per foot. 
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(a) Resistance- thermometer probe. 



45 v d.c . 



Transformer voltage ratio = 1:31.6 


(b) Resistance- thermometer circuit. 


Figure 13. - Re si stance- thermometer instrumentation. 
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(a) Oscillogram 
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(b) Graphical construction. 


Figure 14. - Determination of time constant of resistance- 
thermometer wire. 
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(a) Thermopile probe and reference thermopile. 


Thermopile probe 



NACA - Langley Field, Va. 



